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Abstract 

The spectral problem on a periodic domain with cracks is studied. An 
asymptotic form of dispersion relations is calculated under assumption 
that the opening of the cracks is small. 

Let il C R 2 be a connected smooth domain satisfying the following condi- 
tions: 

• fl is invariant under the shift (x, y) i— ► (x + 1, y), 

• the domains fi„ :=ftn ( [n, n + 1] x R J , n € Z, are bounded, 

• the points A n :— (n, 0) are interior points of Q. 

Assume that the vertical lines L n — (n, R), n <E Z, are non-tangent to the 
boundary of f2 and denote T^ := L n n fi \ B e (A n ), where B e {A n ) is the ball of 
radius e centered at A n , 




Figure 1. Domain f2 



Consider in L 2 (J7) the operator H e which is the Laplacian with the Neumann 
boundary conditions on <9fi and on all the curves T^. We are studying some 
spectral properties of H e as e tends to 0. 

The operator H £ can be viewed as the Hamiltonian of elementary cells J7 ra 
coupled by small windows. For positive e one has a periodic system with a 
band spectrum, while at e = the system is decoupled. There is a number 
of papers concerning the spectral asymptotics for waveguides with windows or 
cracks in various settings, see e.g. [5 H51I7H 9"] . In particular, the paper [5] studied a 
situation which is close to ours, more precisely, the case of the Dirichlet boundary 
conditions on a similar periodic structure with particular symmetries (straight 
strip with symmetric cuts) , and the asymptotics of the lowest dispersion relation 
was calculated. All the papers cited use the method of matching of asymptotic 
expansions jB], which is quite sensible to the geometric properties. We give 
here another proof based on the analysis of operator pencils involving boundary 
integrals and inspired by the recent monograph [1 . This allows one to obtain the 
asymptotics of the spectrum in rather general situations (the assumption made 
above can be relaxed is various directions, we have just chosen a basic situation 
in order to keep simple notations). Integral operator pencils of a similar type 
were used by Pierre Duclos with co-authors as a part of the skeleton method 
for the study of multi-particle problems [2] • 

By periodicity, the spectrum of H e can be studied with the help of the 
Bloch theory. For 9 € [0, 2n] consider the values E n (9;e) for which there ex- 
ists a non-zero function \E' n (x, y; 9, e) satisfying the equation — A^ n (a;, y; 9, e) = 
E n (9,e) 1 ^ n (x,y;9,e), the above boundary conditions, and the Bloch quasiperi- 
odicity condition ^ n (x + l,y;9,e) = e l8 ^ n {x 1 y]9,e)] here 9 € [0, 2-k] is a real 
number called quasimomcntum. By a suitable reordering one can assume that 
the functions 9 — > E n (9,e) are continuous and 27r-periodic; these function are 
called dispersion relations of H 6 . The image of a dispersion relation is called a 
spectral band of H £ , and the spectrum of H £ is the union of the spectral bands. 

We are interested in the asymptotic form of the dispersion relations for 
small e. By the standard variational arguments, see e.g. Proposition 1.3 in [9 
or proposition 9.8 in [3], each function E n (-,e) converges to a constant function 
equal to an eigenvalue of the decoupled system, i.e. to an eigenvalue of the 
Neumann Laplacian on Qq. Our main result is the following theorem. 

Theorem 1 Let E be a simple eigenvalue of the Neumann Laplacian in D,q 
and u be the corresponding normalized eigenfunction. As e —¥ 0, the dispersion 
relation of H e near E is of the form 

E(9,e) =E+ J*?-\u(Ao) - e l9 u(A 1 )\ 2 + o(-^-) 
lose ' V off fJ 



loge|' Vlog 2 e; 



uniformly in 9 € [0, 2tt]. 



We note that E can be any simple eigenvalue, and not just the bottom one. 
The rest of the paper is devoted to the proof of theorem [1] 
Let G n (r,r'; z) be the Green function of the Neumann Laplacian A^„ in f2 n , 

i.e. the integral kernel of the resolvent (N n — z)~ l . Recall that, if z is not an 

eigenvalue of N n , then the boundary value problem 



(-A- 2 )/ = 0inO„, |£ 



en n = ff ' 



where v is the exterior normal vector, has the unique solution 

f(r)= I G n (r,r'-z)g{r')dl{r'). 
Jda„ 

Consider a Bloch solution fy of H e corresponding a quasimomentum 8 and 
to an eigenvalue z = E(9). Introduce the functions 

Uy):= d -^ , ye(-e,s). 

OX x=n 

Denote by *$> n the restriction of Vf to il n - One has, obviously, 

®n(x,y) = - / G n {x,y;n,y';z)f n (y')dy'+ G n (x,y;n+l,y'; z)f n+1 (y')dy' . 

The Bloch condition for ^f implies f n +i{y') = e te f n (y'). The continuity of ^ 
at x — n takes the form ty n (n,y) = 5'„_i(n, y), or, using the above integral 
representation, 



G n (n,y;n + l,y';z)f n (y')dy'- / G n (n,y;n,y'; z)f n (y')dy' 

■ J — e 

G n -i(n,y,n,y';z)f n (y')dy' - e~ lt} \ G n -i(n,y;n- l,y';z)f n (y')dy'. 



Using the obvious identity G n (x,y]x' ,y'; z) = Gq(x — n,y;x' — n,y';z) one 
arrives at a single integral equation 

Kg tS {z)fe,e ■■= J K e , £ (y,y';z)fg^(y')dy' = 0, fe, e ~ fo, (1) 

with the integral kernel 

K e , e (y, y'; z) = G (0, y; 0, y'\ z) + G (l, y; 1, y 1 ; z) 
- e ie G (0, y; 1, y'; z) - e^G^l, y; 0, y'; z). 

In order to determine the dispersion relations of H e one needs to find the val- 
ues of z for which the equation ([1} has non-zero solutions, i.e. the nonlinear 
eigenvalues of Kg jB . 

As noted above, integral equations of this type were studied in details in the 
recent monograph [T], and here we add some details needed for the treatment 
of the periodic problem. First of all, by Lemma 5.6 in [I], for each fixed 6, the 
equation ([T]) has a unique solution in a neighborhood of E if e is small enough. 

Like in pQ introduce the Hilbert space 



X e -=W- IMIjf, :=/ ^J e 2 ~ y 2 \ V {y)\ 2 dy < +x>) , 

and the linear space 

y E :=WeCKe]:^eI £ }. 
By the trace theorems, the functions fg e belong to both X £ and Y 6 . 



We emphasize that the symbol (•, •) will always denote the usual scalar prod- 
uct in L 2 [— e,e], i.e. 



(f>s)'=J_ f(y)g(y)dy. 

Using Lemma 5.1 in [I] and the spectral representation of the Green function 
Go one can decompose Kg iE as follows: 

K ^(z):=~L e + -^- + Re, e (z), (2) 

2n E — z 

where L £ is the integral operator 

L e f(x) = / log | a; - y\ f(y) dy, 
Mg tS is the rank one operator 



M e,ef(x) = ug, £ {x) / U6, e {y) f(y) dy, u e . e (y) := u(0, y) - e l u(l,y), 
and R$ <s (z) is an integral operator 

R0,e(z)f(x)= / Re,e(x,y;z)f(y)dy, 



where the kernel Rg e (x,y;z) is holomorph in (z,9) and of Holder class C 1+a 
with respect to (x,y); here a £ (0, 1) is arbitrary. Note that the functions ug iE 
are in both X e and Y e as well. 

By Lemma 5.2 in [I], the operator L e : X e — > Y" e is a bijection for £ small 
enough, and, by Lemma 5.4 in pQ, one has, uniformly in # G [0,27r] and z in a 
neighborhhod of E, the norm estimate 

\\L^Re.M\\Hx s .x c) ^o(^)- (3) 

here and below L(X 6 ,X e ) is the Banach space of bounded linear operators on 
X £ . We will need the following estimates: 

Proposition 2 As e — > one has, uniformly in 9 , 

(^Ve,^) = ^^ + 0(-i-^). (4) 

loge V | log e | ^ / 

Futhermore, there exists A > «rat/i i/ie following property: If B is a bounded 
operator in X e then, as e — > 0, one has 

|(BL E u 9 , e ,Ufl, e )| < p^j . (5) 

Proof As in Lemma 5.4 of [T] we use an explicit form for the inverse of L e and 
some properties of the finite Hilbert transform. We have 



1 f £ y/e 2 - y 2 u' e e {y) 

L F Ug Ax) = ; / : dy 

E ' V ' ir 2 V^^J-e x-y y 

7T log I Ve 2 - x 2 



where 

a(ue, s ) = ue, e (x) - L e v 9 ^(x) 

(the choice of x is arbitrary) with 

7r 2 Ve 2 - x 2 J-e x — y 

here and below all the integrals are understood in the sence of the Cauchy 
principal value. 

We will use the following well known estimate: for any a € (0, 1) there exists 
C > such that for any ip e C a [— 1, 1] and x € [—1,1] there holds 

r 1 v(p) 



x-y 



dy < C\\<p\\c«[-i,i]- 



One has then 
1 / A 6 

< — 



y/e 2 — x" 1 



y/e 2 - y 2 u' ge {y) 



x-y 



l-E %/e 2 — x 



-dx 



dy 



T 3/2 



Ve 2 - y 2 u' g , e {y) 
x-y 

f 1 V^-y 2 u'e, e {ey) 



\ 2 
dy] dx 



£<»[—£,, 



< 



-l 

Ce 

^3?2 



x-y 



dy 



V 1 -y 2u 'eA £ y) 



£°°[-l,l] 



C«[-l,l]' 



(6) 



and one obtains ||i"i||x e = 0(e) uniformly in 9. This implies (Ii,ug, e ) = 0(e) 
as well. Hence it is sufficient to study the asymptotics of (I2, ug >e ). We have 



it \ 1 / r \ueA*)\\ r 



LeV6,e{y)—j===dy 

v £ - y 



^log|e/2 



(Ki - K 2 ). 



There holds 



\ue,e{x)\ 
y/e 2 — x- 



-.dx = 



' l \ue, e {ex)\ 

-1 yr^ 



-dx 



1 ,mi2 /" da; /" 



\U9,e 



1 2 1 1 

(ex)| - |ue, £ (0)| 



-dx 



yfY^x 2 
-n\u e . e (0)\ 2 +0{e) 1 



where the remainder term is uniform in 6. 

To estimate K 2 let us note first that, by the same arguments as in ([5]), 



f V x-y 2 u ' e - M dV ^ KZ{X) = ° {£) 



uniformly in x and 9. Therefore, 

1 f 6 

sup \L £ vg. £ (x)\ = — sup / log \x - y\ ■ K 3 (y) 

xe[— s,e] K xe[—e,s] J —e 



< — sup \K 3 (y)\ sup 

ft I/e[-e,e] i6[-e,e] 



log|jc-j/| 



\A 2 - 2/ 2 



dy 



r/y 



uniformly in 0. Finally, 



- sup 1^3 (y) I log- = 0(e log e) 
ft ye[-e, e ] l 



I a 2 1 < sup \L s ve,e{x)\ ' ' ay = 0(£l0g£) 

ccef-e,el ./-e \/£ 2 - ?T 



-e V^ - 2/ 2 



We obtain 



(L e l Ug t e,U0 tE ) = (h,Ug. £ ) + (I 2 ,U0 <e ) = (Ii,u g , € ) + 



Ki 



K?. 



| 7r|u e , £ (0)| 2 + O( £ ) | Q( £ log £ ) _ |^, £ (0)| 2 | f> 



Trlogf Trlogf 
1 



7T log ; 



7rlog ■ 



loge 



log e | 2 



which proves ((!}. 

To show ([5]) we use first the Cauchy- Schwartz inequality for X £ , 



(BL £ 1 Ug }£ ,Ug, £ )\ = / \/ E 2 - X 2 ■ BL £ l UB, e {x) 



ue, s (x) 

\/e 2 — x 2 



,dx 



< \\B\\ L{Xe ,x t ) ■ WL^uoAx. 



U0, £ (x) 



e \j£ 2 — X" 



-.dx 



< a\\B\\ L (x e ,x c )\\L £ l u e ,e\\x e , a > 0, 

and the previous estimates show that ||L~ 1 ue !£ ||x e = 0(1/ loge). Proposition 
[2] is proved. 

Let us conclude the proof of theorem [TJ Assume that Kg^(z)fg t e(z) = 
with fg >e y^ and z — E(9, e) sufficiently close to E. Using the invertibility of 
L £ and the decomposition 1(5} one obtains 



(l-2irL- 1 Rg^f e 



2tt 



(u0,e,fe,e)L £ 1 Ug y 



E- z 
The operator 1 — 2irL~ 1 Rg^ is invertible for e small enough, therefore, 

2tt 



fo 



E- z 



(l-2irL £ 1 Rg. £ ) (ug, £ ,fg. e )L e 1 Ug t , 



which implies (ug }E ,fg t£ ) ^ (otherwise one would get fg_ £ = 0). Taking the 
scalar product (•, •) of the both parts with ug £ one arrives at 



= E{9,e) =E-2n((l- 2nL; 1 Re, £ ) L^ug^ug,,). 



Using ([3]), let us fix B > such that 



\\L- 1 RoAz)\\l { x^) < jj^i (7) 



for small e. As the map / n- (f,ug_ e ) is continuous in the A £ -norm, we have, 

2irB 

log 5 



for sufficiently small e such that ,^ B , < 1, 



(L~ l Ue, e , Ug. e ) + ^2 ( ^L^i?^) L~ l U$, e , Ug. 5 \ =". D x + D 2 . 



fe>l 

2^1 



In view of the estimate (@} for Di, it is sufficient show that £>2 = 0(1/ 1 loge 
By © and ©, 

(2nL 7 R e , e ) L 7 u $ , e , U9tE ) < —(— j , 

hence 

which finishes the proof of theorem [TJ 
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